Solitary versus Shock Wave Acceleration in Laser-Plasma Interactions 
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The excitation of nonhnear electrostatic waves, such as shock and solitons, by ultraintense laser 
interaction with overdense plasmas and related ion acceleration are investigated by numerical sim- 
ulations. Stability of solitons and formation of shock waves is strongly dependent on the velocity 
distribution of ions. Monoenergetic components in ion spectra are produced by "pulsed" reflection 
from solitary waves. Possible relevance to recent experiments on "shock acceleration" is discussed. 



In the interaction of superintense (/ ~ 10^® — 
10^1 W cm~^) laser pulses with overdense plasmas (i.e. 
having electron density Uf. > ric, the cut-ofF density), 
the light pressure ranges from Gigabar to Terabar values 
and, like a piston, may sweep out and compress the laser- 
produced plasma pushing its surface at nearly relativistic 
speeds. Such combination of strong compression and ac- 
celeration is often described as the generation of strong 
shock waves (briefly, "shocks") propagating towards the 
bulk of the plasma. These light-pressure driven shocks 
may cause heating of solid targets, as experimentally in- 
ferred [l| and be of interest for the production of warm 
dense matter. In addition, shocks may lead to accel- 
eration of ions if the latter are reflected by the shock 
front as a moving wall. In moderately overdense and hot 
plasmas, where the shock waves are of collisionless na- 
ture, shock acceleration may lead to higher ion energies 
than the widely studied target normal sheath accelera- 
tion (TNSA) mechanism 0] , as suggested on the basis of 
numerical simulations 0] . Recently, two experiments em- 
ploying hydrogen gas jet targets and CO2 lasers reported 
on monoenergetic proton beams produced via shock ac- 
celeration 0, The study of laser-driven shocks in 
the laboratory may also help the understanding of simi- 
lar processes in astrophysical environments [1, and refer- 
ences therein] 

The above mentioned simple picture of ion accelera- 
tion indicates that, as far as the shock velocity Vs is con- 
stant, the reflected ions should have velocity 2vs and pro- 
duce a monoenergetic peak in the spectrum. Silva et al. 
3 report that such peak would evolve into a spectral 
plateau due to further acceleration in the sheath field at 
the rear side of the target, so that suppression of such 
field (for instance by producing a smooth density gra- 
dient at the rear side) could preserve monoenergeticity. 
However, in the present Letter we show that obtaining 
monoenergetic spectra is not straightforward, even ne- 
glecting the side effect of the sheath field. We observe 
monoenergetic peaks in the simulations only when short- 
duration ion bunches are accelerated by solitary waves 
generated by the laser-plasma interaction. Moreover, the 
formation and the evolution of both soliton- or shock-like 
waves is highly dependent on the velocity distribution of 



background ions. A crucial point is that, for a ion dis- 
tribution without a velocity spread, either all the ions 
or none of them may be reflected from a moving electro- 
static field front, depending on the ratio ^^ax/^^s where 
<I>„^x is the potential jump at the front. In the first case, 
a strong loading of the moving structure occurs, causing 
deceleration and collapse of the shock or soliton and a 
consequent broadening of the reflected ions spectrum. 

Before describing our simulation results we briefly re- 
call some key aspects of the fluid theory of collisionless 
electrostatic shocks and solitons The most relevant 
point for what follows is that the existence of reflected 
ions is integral to collisionless shocks and not a secondary 
effect; it may be used as a signature of "true" shock for- 
mation in the simulations. When reflected ions are ab- 
sent, solutions corresponding to electrostatic solitons are 
found. Hence, a necessary and general condition for such 
solitons to exist with a velocity Vg is that the electrostatic 
potential energy jump e$ has a peak value 
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so that background ions are not reflected by the soliton. 
Within the fluid theory with the electrons in an isother- 
mal Boltzmann equilibrium Q at the temperature Te, 
the condition on the potential poses an upper limit on 
the Mach number M = Vg/cs < 1.6, where Cs — ^/T^jmi 
is the speed of sound. The other condition is that the 
soliton must be supersonic, i.e. M > 1. 

Since we focus on basic aspects we restrict to one- 
dimensional Particle-In-Cell simulations for the sake of 
simplicity and high numerical resolution. We found that 
the numerical results converged only for sufficiently high 
values of the number of particles per cell N^, although 
qualitatively similar features were observed also for lower 
A'p. This suggests that kinetic effects play an essential 
role and thus low-density tails in the distribution func- 
tions must be resolved accurately. In the simulations 
reported below, = 800 and the spatial and tempo- 
ral resolution Aa; = cAt = A/400, where A is the laser 
wavelength. 

We now analyze a representative simulation with the 
following set-up. The laser pulse was linearly polar- 
ized with a peak amplitude ao = 16, where = 
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FIG. 1. (Color online). Snapshots of the evolution of a Solitary Acoustic Wave (SAW) at eight different times. The upper, 
middle and bottom rows show the ion density m, the electrostatic field Ex and the contours of the fi{x,px) ion phase space 
distribution in a log^Q-scale. The laser pulse impinges from the left reaching the plasma boundary {x = OA) at t = OT. 
Simulation parameters are ao = 16, r = 4T and Ue — 20nc. rii is normalized to Uc, Ex to Eo — nieujc/e and px to mpC. 



0.85(/A2/101^ W cm"V^l2)^/^ and duration r = 4T 
(FWHM), with T the laser period; the temporal pro- 
file is composed by IT long, sin^-like rising and falling 
ramps and a 3T plateau. The plasma had a slab, square- 
like profile with initial ion and electron densities rii = 
Uf, = riQ = 20nc and 15A thickness. Ions were protons 
[Z = A = 1). For reference, the laser pulse front reaches 
the front surface placed at a; = at the time i = 0. 

Eight snapshots of the profiles of ion density and 
longitudinal electrostatic field E^. and of the fi{x,px) 
phase space distribution of ions are shown in Fig[T] In the 
early stage of the interaction, the laser pulse accelerates 
a fraction of strongly relativistic electrons with energy of 
several nieC^ which penetrate into the target and later 
recirculate across it, driving heating of bulk electrons. A 
solitary structure is generated at the front surface under 
the action of the laser pulse and then propagates into the 
plasma bulk at a constant velocity Vg ~ 0.05c. This value 
is close to the "hole boring" velocity v^.^, [1, |^ which in 
our case is given by 

y - a c f ^— i-^ — — — V^^ (2) 
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where the reflection coefficient is measured to be i? ~ 
0.75 in the simulation of FiglU yielding — 0.06c. At 
t = 65r, the solitary structure is located at x ~ 3.7A 
(first frame of Fig[T]). The ion density has a very strong 
spike, reaching values up to ~ 9 times the background 
density. The electric field around the density spike has a 
sawtooth shape. For reference we call the structure we 
observe a SAW (for Solitary Acoustic Wave). The snap- 
shots of SAW density and field profiles are qualitatively 
resemblant to those of electrostatic solitons as described 
in Ref.Q. However, several of the dynamics features we 
describe below are not described by standard fluid theory. 

Numerical integration of over x for the SAW at t = 
65T (first snapshot in FigH]) yields e$,„,, ~ 0.78i;oA ~ 



4.9meC^. Since ■miv'^/2 ~ 2.3meC^ < e^,^^^ one would 
expect ion refiection to occur promptly from the SAW 
field. This occurs indeed in the simulation, producing 
a bunch of high-energy ions with very small momentum 
spread as can be observed at t = 70T and t = 75r. 
However, violation of the condition ([T|) and consequent 
ion reflection do not destroy the SAW; at t ^ 75T, the 
SAW is still present with almost unperturbed velocity, 
and ion reflection has stopped. Between t = 75T and 
t = 85T a second bunch is formed in a very similar way, 
as observed between the third and fifth snapshot in FiglT] 
Both bunches correspond to a monocncrgctic high-energy 
peak in the spectrum, as shown in Figl2]a); the FWHM 
of the peak at t = 75T is less than 1%. 

The acceleration of a third bunch at t ^ 90T is cor- 
related with a much stronger perturbation of the SAW 
whose amplitude and velocity decrease. As a conse- 
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FIG. 2. (Color online), a): ion spectra from the simulation of 
FigH] at the times t = 75T, 90T and lOOT. Only ions in the 
same spatial range as FiglU i-e- 3-5 < a; < 7.5, are included 
in the spectra, b): temporal evolution of the maximum (red 
line) and the minimum value (blue line) of the electric field 
for the SAW structure observed in Fig[l] The vertical dashed 
lines mark the instants t — 65T, 75T and 90r at which ion 
reflection events start as seen in Fig[T] 
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FIG. 3. (Color online). Snapshots at three different times 
of the Wi, Ej; and fi{x,px) for a simulation identical to that 
of Fig[T]but for an initial ion temperature Ti ~ 10~^meC^ = 
5 keV, showing the onset of "continuous" ion reflection. The 
scales are the same as in Fig[T]to show that the perturbation 
in Hi and Ex has a lower value in the present case. The insets 
show the corresponding ion spectra, limited to the ions in the 
region around the shock and excluding ions accelerated at the 
target boundaries. 



quence, the high-energy part of the spectrum broadens 
up, and monoenergetic features are lost as seen in Figl2]a) 
at t = lOOT. 

The apparent behavior of the SAW is related to the 
observation that the electric field amplitude is not con- 
stant in time, as for the "standard" soliton solution, but 
it oscillates as shown in Figl2l The temporal behav- 
ior of the maximum and minimum values shows that 
the electron cloud around the ion density spike oscil- 
lates back and forth. The oscillation is quenched after 
the generation of the second fast bunch at t ~ 75T. 
At this instant, the potential jump at the SAW front is 
e$„-,ax — 0.5EqX ~ S.lmgC^, slightly above the stability 
threshold ([T]). Eventually, the overall amplitude strongly 
decreases after the generation of the third bunch, as also 
observed in the last frames in FigH] Apparently, such 
"collapse" of the SAW occurs after a "collision" with a 
slower, counterpropagating structure that can also be no- 
ticed in the frames of Fig[T] 

The above reported observations led us to infer that 
the local velocity distribution of ions has a crucial role 
in determining both the stability of solitary waves and 
the possibility to generate "true" shocks whose signature 
would be, in the framework of fluid theory, the forma- 
tion of a continuous flow of reflected ions. If the ions 
are "cold" , i.e. have no energy spread, reflection from a 
moving potential barrier may occur cither for none or for 
all of the ions. In the latter case, the wave may quickly 
lose its energy by accelerating the whole bulk of ions. 

In order for a collisionless shock to form, it must be 
possible for the wave to "pick up" from the ion distri- 
bution only a fraction of the ions with the right energy 



to be reflected. This can occur in a plasma with energy 
spread, i.e. warm ions. Fig|3] shows results of a simu- 
lation with identical parameters as Figdl but the initial 
ion temperature = 5 kcV. Formation of a shock-like 
structure with steady reflection of ions is observed. The 
shock front velocity is found to decrease in time which 
may be interpreted as due to the wave energy transfer 
into reflected ions. As a consequence, with respect to 
the "cold" ions case narrow monoenergetic peaks almost 
never appear at any time in the spectrum (see the insets 
in Fig|3]) and the high-energy tail progressively broadens 
towards lower energies, although the cut-off energy is a 
bit higher with respect to the "cold" case. Thus, in this 
spectral "plateau" is produced because of shock 
deceleration, rather than due to further acceleration of 
the reflected ions in the rear side sheath. Some modula- 
tions are apparent both in the phase space density and 
in the spectral plateau of reflected ions; they may be re- 
lated to the oscillation of the electric field at the wave 
front, similar to FiglJJ 

The effect of the background ion distribution can be 
also evidenced studying a case with parameters identical 
to Figd] but a shorter length of the plasma (3A). In 
this case, not shown here for briefity, the SAW generated 
at the front side of the plasma eventually reaches the 
rear side and propagates in the expanding sheath until it 
reaches the region in which the ion velocity is such that 
the ions are now reflected by the SAW potential (since 
the ions in the sheath move in the same direction of the 
SAW, their velocity in the rest frame of the SAW is lower 
than in the "laboratory" frame, hence they are reflected 
more easily). From this point on the SAW quickly loses 
its energy in accelerating ions and collapses. A similar 
mechanism was also described by Zhidkov et al. [lol | . 

When the laser pulse duration is increased with respect 
to the above reported simulations, we observe multiple 
peaks of the ion density and sawtooth spatial oscillations 
of the electric fleld. Although a structure showing multi- 
ple oscillations behind a front is reminiscent of a collision- 
less shock wave as seen in textbooks, again we do not ob- 
serve in general a steady ion reflection at the front. The 
density peaks move at different velocities and thus dis- 
perse in time. Hence the structure may be interpreted as 
a multi-peak SAW, generated due to pulsed hole boring 
acceleration Q at a rate which is approximately the same 
for simulations having same pulse intensity and plasma 
density, so that a longer pulse duration allows accelera- 
tion of a sequence of ion bunches. 

A representative "long pulse" simulation is shown in 
FigHl The laser pulse has peak amplitude oq = 16 (like 
the above reported "short pulse" simulations) with 5T 
rise and fall ramps and a GOT plateau. The plasma den- 
sity and thickness are lOric and 15A, respectively. These 
parameters are close to those of previously reported ID 
simulations (see Fig.l in Ref.Q) and indeed some fea- 
tures observed in the ion phase space and density profiles 
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FIG. 4. (Color online). Snapshots of multi-peak structures 
at three different times from a "long pulse" simulation. The 
inset shows the spectrum at t = 45r, limited to the ions in 
the same region as the plots. The vertical dashed line marks 
the initial position of the rear side of the target. Simulation 
parameters are ao = 16, t — 65T and Ue — lOric. Units are 
the same as in Fig[T] 



look very similar. However, in our simulations the higher 
mimber of computational particles allows us to highlight 
additional details in the phase space distribution of Fig|4] 
(bottom row) such as vortex structures behind the front, 
corresponding to "trapped" ions bouncing between ad- 
jacent peaks, where a potential well is formed. The ob- 
servation of ions trapped in the multi-peak, nonlinear 
structure of the electric field and the quite broad distri- 
bution along the momentum axis of the most energetic 
ions suggest that those latter may be actually acceler- 
ated by "surfing" the longitudinal wave structure, rather 
than being merely reflected by the field at the front. In 
addition, we observe significant oscillations of the electric 
field ahead of the wave front (i.e. of the rightmost density 
peak in Fig|4]), which may have been excited by fast elec- 
trons and are related to modulations in the momenta of 
the highest energy ions. The spectra of the latter is quite 
broad (inset of FigU] even neglecting the contribution of 
TNSA at the rear side. 

Our study suggests that generating highly monoener- 
getic ions by "shock acceleration" is not straightforward, 
as in our simulation we observe narrow spectra only as 
resulting from SAW "pulsations" and as a transient ef- 
fect, since as shown in FigsH] and [2] the SAW collapse 
ultimately produces a broad spectrum masking the mo- 
noenergetic peak. Producing a monoenergetic spectrum 
might not coexist with efficiency, because the reflection of 
a large fraction of ions by the moving structure (either a 
shock or a soliton) would ultimately cause a strong load- 
ing effect decreasing its field and velocity. This is in qual- 
itative agreement with the experimental results of Haber- 
berger et al. where the number of ions (^ 2.5 x 10^) in 
the narrow spectral peak at 22 MeV implies a conver- 
sion efficiency < 10^^ of the 60 J pulse energy. A direct 



comparison of our simulations with the experiment in 
Ref.[5|, which involves several additional issues such as 
e.g. the pulse shape and the plasma profile, is beyond 
the scope of the present Letter. 

All of the above described effects related to SAW gener- 
ation, dynamics and acceleration of ions in the bulk are 
observed only with linear polarization, for which elec- 
trons are heated up to relativistic energies. For circular 
polarization, no ion acceleration occurs in the bulk. Al- 
though Palmer et al. 01 report on "protons accelerated 
by a radiation pressure driven shock" , we find that the 
experiment can be better described as evidence of hole 
boring (HB) acceleration [§, Nevertheless, the two 
different mechanisms might produce similar ion spectra 
because the velocity of the shock launched at the sur- 
face is of the order of f ht [Eq. ^] , while the hole boring 
process produces via wave breaking ion bunches of veloc- 
ity Vi ~ 2vhb 01 if ions were "reflected" from the 
plasma surface as a piston driven by light pressure. 

In conclusion, in our simulations we observe the 
laser-driven generation of nonlinear, collisionless acoustic 
waves having either a solitary or a multipeak structure 
depending on the laser pulse duration. In particular we 
observe a "pulsation" of solitary waves leading to the 
generation of small numbers of monoenergetic ions. In a 
cold ion background, wave loading effects prevent "true" 
shock wave formation and efficient monoenergetic accel- 
eration. In general the dynamics of "shock acceleration" 
in the plasma bulk appears to be more complex than the 
simple picture of "reflection from a moving wall" . 
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